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Abstract. We discuss two approaches to searches for gravitational-wave (GW) and electromag- 
q ' netic (EM) counterparts of binary neutron star mergers. The first approach relies on triggering 
archival searches of GW detector data based on detections of EM transients. We introduce a 
quantitative approach to evaluate the improvement to GW detector reach due to the extra 
information gained from the EM transient and the increased confidence in the presence of a 
m : signal from a binary merger. We also advocate utilizing other transients in addition to short 
gamma ray bursts. The second approach involves following up GW candidates with targeted EM 
[•j - ] ' observations. We argue for the use of slower but optimal parameter-estimation techniques to lo- 
|_ , calize the source on the sky, and for a more sophisticated use of astrophysical prior information, 
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Adv anced groun d- based gravitational - wave detectors LIGO Harry fc the LIGO Scientific Collaboration 
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including galaxy catalogs, to find preferred followup locations. 
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1. Introduction 



2010I) ] and Virgo [Virgo Collaboration! (2009)] are expected to begin taking data around 



t-H ■ 2015, operating at a sensitivity approximately a factor of 10 better than their initial 
counterparts. Gravitational waves (GWs) emitted during the late inspirals and mergers 
of compact-object binaries will be one of the main sources for these detectors. Although 
predictions for merger rates of bina ry neutron stars and neutron-sta r - black-hole bi- 
naries are highly uncertain (see, e.g, Mandel fc O'Shaughnessv ( 2010l) ). we may expect 



detections at a rate between once per few years and a few hundred per yea r, with perhaps 



a few tens of detections per year being most likely [Abadie et al.l (|2010al) ]. 

These detections would usher in an era of genuine gravitational-wave astronomy, with 
GWs being used as another tool to observe the sky. There has been a lot of discussion 



of the promise of multimessenger observations in the literature (see, e.g., iBloom et al 



(2009)). [n fact, the recent focus on multimessenger astronomy has been so exclusive 



that it is worthwhile to remind ourselves that a significant amount of astrophysics can 
be extracted from the GW observations alone, since the GW signal encodes the masses 
and spins of the binary components, and can be used to probe astrophysics, strong- 
hold gravity, and cosmology even in the absence of electromagnetic (EM) observations 
of counterparts to GW events. Nonetheless, there is no doubt that observing both EM 
and GW counterparts of the same event would be of great astrophysical significance, and 
would allow us to settle crucial questions such as the origins of short GRBs. 

In general, such observations can be achieved in three different ways: (i) through en- 
tirely independent, serendipitous observations of GW and EM signatures of the same 
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event, (ii) through triggered searches of archival GW data based on EM transients ob- 
served during surveys, or (iii) through telescope pointing to search for electromagnetic 
followups of GW candidates. Here, we describe some thoughts about future possibilities 
for triggered GW searches based on observed EM transients, and possible improvements 
to recently started efforts to follow up GW triggers with targeted EM observations. 

2. GW searches triggered on EM transients 

The detection of an electromagnetic transient which may originate from a compact- 
object binary merger will increase the a priori probability that a given stretch of data from 
the LIGO- Virgo ground-based gravitational- wave detector network contains a signal from 
a binary coalescence. Additional information contained in the electromagnetic signal, such 
as the sky location or distance to the source, can further help to rule out false alarms, 
and thus lower the necessary threshold for a detection. 

The LIGO Scientific Collaboration and the Virgo Collaboration have previously used 
short, hard Gamma Ray Bursts ( SGRBs) as tr i ggers f or searches for compact binary 
coalescences in GW detector data jAbadie et al. ( 2010bh ] . SGRBs are believed to be as 



sociat ed with relativistic jets from BNS or NS-BH mergers (see, e.g. jLee fc Ramirez-Ruiz 



(|2007l) ). However, the impact of the SGRB observation on gaining confidence in the pres- 
ence of a GW signature from a binary merger has not been quantitatively estimated. In 
fact, most of the observed SGRBs originate too far away to detect the associated GW 
signal: the average luminosity distance for the 16 SGRBs wi th confident host identifi- 
cations and redshift measurements compiled by iBereerl (l2010h is approximately 5 Gpc, 



which exceeds the anticipated event horizon of Advanced LIGO by more than a factor 
of 10. Hence, even when an SGRB is detected, it is unlikely that a gravitational-wave 
counterpart is observable, so the detectability threshold is not lowered by as much as 
c ould be anticipated. 



Kellev et al.1 (|201lh develop a Bayesian framework to accurately incorporate the ad- 
ditional prior information on the existence of a signal and its timing, sky location, and 
potentially distance and inclination, in order to quantify the benefits of triggered searches 
in improving the detector reach. The odds ratio for the GW detector data d to contain 
a signal relative to the noise-only model is 

p(GW|d) _ p(GW)p(d|GW) 
p(N|d) " " p(N) p(d|N) ' 

where the first term is the a priori probability of having a GW signal present, and the 
second term is known as the Bayes factor, 

B = [ P (6\GW) ■ e <d\h{e)>-\<h{e)\h(9)> d ^ (2 2) 



where the evidence integral is computed over all possible values of the parameters 6 of the 
signal model with their corresponding priors p(0|GW). We can approximate the Bayes 
factor by assuming that the likelihood in the previous integral is roughly constant over 
a fraction r\ of the parameter space for which the signal model is supported by the data, 
where it has a simple dependence on the signal-to-noise ratio (SNR), and zero elsewhere: 

»«f (2.3) 

The first term in Eq. (|2.3[) can be increased by the detection of an EM transient, 
depending on our a priori confidence that this kind of transient is associated with a binary 
coalescence, the probability that the coalescence occurs within the detection volume, and 
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the accuracy of the merger time reconstruction from the transient. The second term is 
also affected by the additional information from the transient on the sky location of the 
source, as well as possibly the distance to the source and its inclination, which constrains 
the priors on individual parameters and increase the fraction of the parameter space 
which offers support for the signal model. Assuming a fixed threshold on the odds ratio 
required for detection (corresponding to a fixed false alarm probability) , the EM transient 
can therefore reduce the SNR threshold for detection by the following fraction: 



CsNR 



SNR EW 
SNR 



In O 



P (GW|EM) 
p(N|EM) ' 



In O 



p(GW) n 



(2.4) 



Preliminary results Kellev et al. ( 201lh ] suggest that optical signatures from r- processes 
in the tidal tails (also known as "kilonovae") could be one of the most promising triggers. 
They are likely to be sufficiently numerous in all-sky surveys like the LSST, differentiable 
from other transients based on the light curve profiles if the observational cadence is suf- 
ficiently high, and could yield moderate improvements in the detectability thresholds. 

The above formalism assumes that optimal coherent searches are used both in the 
presence and in the absence of an EM trigger. If the untriggered search is a subopti- 
mal inco herent search while the tr iggered search is optimized because of the known sky 
location [Harry fc Fairhursti f)201 lh ] . the improvement can be even more significant. 



3. EM followups of GW candidate events 

The benefits o f following up GW candidates with targeted EM observations have long 
been recognized Finn et al. ( 1999h ]. Re cently, the first pilot program for following up 
GW candidate events has been activated Abadie et al.l (|201 lh ] . There are two significant 
complications in this effort. One is the speed with which GW data can be analyzed 
to yield information about detection confidence and the sky location of GW candidates. 
The other is the large positional uncertainty associated with GW detections, which could 
encompass tens or even a few hundred square degrees on the s ky depending on the S NR 
of the candidate and the detector network configuration (e.g., Nissanke et al. ( 201lh ). 

A significant amount of effort has been expended to allow for very rapid processing of 
GW data with the aim of achiev ing latencies o f only a few seconds or tens of seconds 
to identify GW candidates (e.g., Cannon et al. ( 201l[ )V Timing triangulation between 
different GW interferometers co mprising the net work has been used so far to rapidly 
localize the source on the sky fe.g. jFairhurstl ( 2009I )). Such incoherent use of detector data 
allows for rapid analysis, but is likely to yield suboptimal results. The "need for speed" 
can be over-stated: it is simply not necessary for many of the possible EM counterparts. 
For example, off-axis optical aft erglows outside the jet opening a ngle will only peak on 
time scales greater than ~ 1 day van Eerten fc MacFadvenl (|201 lh ] . We may hope that if 
we are inside the jet opening angle, the GRB will be picked up by a mission like Swift; but 
even if not, the optical signal has decay timescales of at least several hours. Timescales 
for lower-frequency followups, such as radio, are even longer. We thus have the luxury 
of applying slower paramete r-estimation techniques tha t perform a coherent analysis of 
all of the detector data (e.g. Jvan der Sluvs et al. (2008)), a nd should be a ble to improve 
sky localization and better quantify the confidence regions [Avlott et all (|2012h ]. 

Even with slower, optimal data-analysis techniques, the confidence regions are much 
larger than the fields of view of typical followup instr uments. Efforts to reduc e the GW 
sky error box have relied on using a galaxy catalog Nuttall fc Sutton ( 2010h ] . In par- 
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ticular, pixels within the error box have been re - weigh ted by the blue-light luminosity 



of the galaxies contained in them [Abadie et al. ( 201 lh ] - While a reasonable first cut 



this approach could be problematic for three reasons. Firstly, blue-light luminosity is a 
proxy for the star formation rate, but mergers could be significantly delayed relative to 
star formation episodes, and red elliptical galaxies are tho ught to contribute significantly 



to present-day merger rates [O'Shaughnessy et al.l ( 2010h ] . Secondly, galaxy catalogs are 
unlikely to be complete to ~ 400 Mpc (the horizon distance of Advanced LIGO for op- 
timally located and oriented coalescing neutron-star binaries). Finally, some fraction of 
mergers may happen outside galaxies altogethe r, as the progenito r binary could have 
be en kicked out of t he ho sts by supernovae kicks [Kellev et al.f(j2010h ]. An ongoing study 
by IVousden et al~ ( 2012 ) aims to account for these shortcomings by employing more 
sophisticated astrophysical priors. 

Despite the improvements mentioned above, finding an electromagn etic counterpart of 
a GW candidate will remain extremely challenging, as discussed by iMetz ger&Berger 
(| 2 1 lh . Coordinated observing among several facilities may be required to cover the large 
uncertainty region with smaller field-of-view instruments. Another alternative that may 
be worth investigating is the deployment of a network of inexpensive robotic telescopes 
specifically with the goal of following up GW candidates, though it will be difficult 
to detect any but the closest afterglows. While it is impossible to guarantee that an 
EM counterpart to given candidate would be found, we can still strive to maximize the 
probability of a successful follow up, with a view to ensuring that a sufficient fraction of 
triggers are followed up to make at least some multimessenger observations. 
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